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ABSTRACT 

We study the significance of major mergers in driving star formation in the early 
Universe, by quantifying the contribution of this process to the total star formation 
budget in 80 massive (M* > 10 1D M©) galaxies at z ~ 2. Employing visually-classified 
morphologies from rest-frame V'-band HST imaging, we find that 55 ±14 % of the star 
formation budget is hosted by non-interacting late-types, with 27 ±s % in major merg- 
ers and 18 ±6 % in spheroids. Given that a system undergoing a major merger continues 
to experience star formation driven by other processes at this epoch (e.g. cold accre- 
tion, minor mergers), ~27% is a iikely upper limit for the major-merger contribution to 
star formation activity at this epoch. The ratio of the average specific star formation 
rate in major mergers to that in the non-interacting late-types is ~2.2:1, suggest- 
ing that the typical enhancement of star formation due to major merging is modest 
and that just under half the star formation in systems experiencing major mergers is 
unrelated to the merger itself. Taking this into account, we estimate that the actual 
major-merger contribution to the star formation budget may be as low as ~15%. While 
our study does not preclude a major-merger-dominated era in the very early Universe, 
if the major-merger contribution to star formation does not evolve significantly into 
larger look-back times, then this process has a relatively insignificant role in driving 
stellar mass assembly over cosmic time. 
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1 INTRODUCTION 


Understanding the formation of massive galaxies is a central 
topic in observational cosmology. The observed peak in the 
cosmic star formation rate (SFR) at z te 2 (e.g. Madau et al. 
1998; Hopkins & Beacom 2006; Le Borgne et al. 2009) indi- 
cates that a significant fraction of the stellar mass in today’s 
massive galaxies is likely to have formed around this epoch. 
However, the principal mechanisms that created this stellar 
mass remain unclear. Was the star formation driven by vig- 
orous, major-merger (mass ratios > 1:3) induced starbursts? 
Or were processes other than major mergers - e.g. cold ac- 
cretion, minor mergers, etc. - responsible for creating the 
bulk of the stars in today’s massive galaxies, as suggested 


by recent theoretical work (e.g. Keres et al. 2009; Dekel et al. 
2009)? 

Modern surveys that access large UV/optically-selected 
samples of galaxies at z > 1.6 have facilitated the empiri- 
cal study of star formation around z ~ 2 (e.g. Daddi et al. 
2004; Erb et al. 2006; Reddy et al. 2006; Daddi et al. 2007; 
Santini et al. 2009; Hathi et al. 2010; Wuyts et al. 2011). 
Star-forming galaxies at this epoch lie on a star-formation 
‘main sequence’ (e.g. Daddi et al. 2007; Reddy et al. 2012), 
where galaxy SFRs are proportional to their stellar masses 
with a slope of unity (relatively passive galaxies he below 
this sequence). The growing body of observational work on 
these galaxies increasingly suggests that much of the cosmic 
Btar formation at this epoch may be unrelated to the major- 
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Figure 1. Example J+H composite images for the morphological 
classes into which we split our galaxy sample. Galaxies are classi- 
fied into spheroids (top row), major mergers (disturbed systems 
with multiple nuclei and clear, extended tidal features, middle 
row) and non-interacting late- types (bottom row). The images 
are ~45 kpc on a side. Note that the images usually appear bet- 
ter on screen than in print. 


merger process. Integral-field spectroscopy of star-forming 
galaxies around z ~ 2 has revealed a high fraction of sys- 
tems with properties indicative of turbulent disks and only 
a modest incidence of major mergers (e.g. Forster Schreiber 
et al. 2006; Genzel et al. 2008; Shapiro et al. 2008; Forster 
Schreiber et al. 2009; Cresci et al. 2009; Genzel et al. 2011; 
Mancini et al. 2011, see also Law et al. 2009, van Dokkum 
et al. 2011). Imaging studies, that have explored the rest- 
frame UV and optical morphologies of star-forming galaxies 
at these epochs (e.g. Lotz et al. 2006; Forster Schreiber et al. 
2011; Law et al. 2012a), have also indicated a preponderance 
of non-merging systems amongst high-redshift star formers, 
suggesting that the role of major mergers may indeed be 
subordinate to that of other processes (such as cold flows or 
minor mergers) in driving star formation in massive galaxies 
at this epoch. 

In a recent study, Rodighiero et al. (2011) have shown 
that ‘starbursts’ - systems that show enhanced star formar 
tion and lie off the main sequence of normal star-forming 
galaxies - have a relatively minor role at this epoch, account- 
ing for around 10% of the cosmic star formation activity. 
However, starbursts can be driven either via major merg- 
ers or t»y dense nuclear star-forming regions (e.g. Di Matteo 
et al. 2007; Dekel et al. 2009; Daddi et al. 2010). More impor- 
tantly, many major mergers share the same star-formation 
characteristics as normal star-forming galaxies at this epoch 
(e.g. Law et al. 2012a, see also Di Matteo et al. 2007, Kaviraj 
et al. 2012) and thus lie on the main star-formation sequence 
itself. As a result, a unique one-to-one mapping is unlikely 
to exist between major mergers and starbursts. To probe 
the relative significance of major-merger-driven star forma- 


tion at z ~ 2, it is desirable to quantify the proportion of 
the total star formation budget that is attributable to sys- 
tems that are morphologically selected as major mergers at 
this epoch. This has not been directly addressed by previous 
work and represents both a quantitative empirical result and 
a useful constraint on theoretical models at high redshift. 

Deep near-infrared imaging from current WFC3 sur- 
veys - which trace rest-frame optical wavelengths at z ~ 2 
— enables us to morphologically classify massive galaxies at 
this epoch and study how star formation activity is appor- 
tioned in terms of galaxy morphology (e.g. major mergers, 
non-interacting late-types, etc). It is worth noting, however, 
that a system undergoing a major merger at z ~ 2 continues 
to experience star formation driven by gas inflow via other 
processes such as cold accretion and minor mergers (major 
mergers can be thought of as simply the 'clumpiest’ part of 
the material flowing in along the cosmic web). Simulations 
indicate that star formation due to these other processes is 
significant at this epoch and possibly comparable to major- 
merger-driven activity (e.g. Dekel et al. 2009). This appears 
consistent with recent empirical work (e.g. Kaviraj et al. 
2012; Law et al. 2012a) which indicates that SFRs in non- 
interacting systems can be similar to those in major mergers. 
Hence, in addition to splitting the star formation budget by 
morphology, it is necessary to consider the fraction of star 
formation in major-merging systems that is unrelated to the 
merger itself (and subtract this from the star formation frac- 
tion hosted by systems with major-merger morphology). 1 

Here, we probe these questions using a complete, rest- 
frame optically-selected sample of massive (M, > 10 10 M 0 ) 
galaxies at z ~ 2, drawn from the WFC3 Early Release 
Science (ERS) programme, which provides unprecedentedly 
deep near-infrared HST imaging and ten-filter photometry 
in the GOODS-South field. Section 2 describes the galaxy 
sample that underpins this study. In Section 3, we describe 
the derivation of galaxy properties e.g. SFRs, stellar masses 
and internal extinctions. We study the proportional contri- 
bution of major mergers to the total star formation bud- 
get in Section 4 and summarise our findings in Section 5. 
Throughout, we use the WMAP7 cosmological parameters 
(Komatsu et al. 2011) and present photometry in the AB 
magnitude system Oke & Gunn (1983). 

2 GALAXY SAMPLE AND 

MORPHOLOGICAL CLASSIFICATIONS 

The WFC3 ERS programme has imaged ~45 arcmin 2 of the 
GOODS-South field in the WFC3 UVIS (F225W, F275W, 
F336W) and IR (F098M [Y], F125W [J], F160W [H]) chan- 
nels, with exposure times of 1-2 orbits per filter. The ob- 
servations, data reduction, and instrument performance are 
described in detail in Windhorst et al. (2011). Together with 
the existing ACS BViz imaging (Giavalisco et al. 2004), the 
data provide 10-band panchromatic coverage over 0.2 - 1.7 
fim, with 5 <7 point source depths of AB ^ 26.1 — 26.4 mag 
in the UV and AB £ 27.2 — 27.5 mag in the IR. 

Here, we focus on an //-band selected sample of 80 ERS 


Note that the situation is significantly different at low redshift, 
where gas-rich major mergers can enhance star formation by or- 
ders of magnitude (e.g. Mihos & Hemquist 1996), because secular 
processes drive star formation weakly. Almost all the star formar 
tion in low-redshift major mergers is, therefore, attributable to 
the merger itself. 


© 0000 RAS, MNRAS 000, 000-000 



The insignificance of major mergers at high redshift 3 


galaxies, that have stellar masses M» > 1O 1O M0 and pho- 
tometric redshifts (calculated using the EAZY code, Bram- 
mer et al. 2008) in the range 1.9 < z < 2.1. The accuracy 
of the photometric redshifts at this epoch is A z ~ 0.1 and 
the nominal time interval defined by this redshift range is 
~0.3 Gyr. The ff-band traces rest-frame V at z ~ 2 and 
the galaxy sample is complete within these stellar mass and 
redshift ranges (Windhorst et al. 2011). Thus, we are not bi- 
ased against galaxies with low star formation rates. Study- 
ing the massive end of the galaxy population restricts us 
to systems that are both bright ( H(AB ) < 24.2 mag) and 
extended which facilitates reliable morphological classifica- 
tion. The narrow redshift interval minimises morphological 
K-corrections and overlap between the spheroid and major- 
merger morphological classes (as we discuss in Section 4). 

Here we classify galaxies via visual inspection of their 
composite J+H images. Since the J and H filters correspond 
to the rest-frame optical wavelengths at z ~ 2, these images 
trace the underlying stellar populations in each galaxy and 
not just the UV-emitting star-forming regions. Visual clas- 
sification of morphologies in the high-redshift Universe has 
been commonly employed in the literature, using rest-frame 
optical HST images that have similar or fainter surface- 
brightness limits compared to the ERS images used here 
(e.g. Windhorst et al. 2002; Cassata et al. 2010; Kaviraj 
et al. 2011; Cameron et al. 2011; Kocevski et al. 2012; Law 
et al. 2012a). Visual classification offers better precision and 
consistency than morphological parameters (such as CAS, 
M 23 , Gini coefficient, see e.g. Abraham et al. 1996; Conselice 
et al. 2003; Lotz et al. 2004; Taylor-Mager et al. 2007), which 
car. be more sensitive to image resolution and signal-to-noise 
(e.g. Lisker 2008; Kartaltepe et al. 2010) but are valuable 
for classifying large datasets, where visual classification is 
prohibitively time-consuming. 

Galaxies are classified into the following three broad 
morphological classes: [1] spheroids [2] non-interacting late- 
types and [3] major mergers, which are disturbed systems 
that exhibit multiple nuclei and clear, extended tidal fea- 
tures. The number fractions in classes [1], [2] and [3] are 
19%, 39% and 42% respectively. Figure 1 presents examples 
of objects drawn from each morphological class. 

3 PARAMETER ESTIMATION: STELLAR 

MASSES AND INTRINSIC STAR 

FORMATION RATES 

A variety of methods have been employed in the literature 
to derive galaxy SFRs. Calibrations (often based on samples 
of local galaxies) can be used to convert X-ray, UV, infrared 
or radio luminosities into estimates of SFR (see e.g. Reddy 
& Steidel 2004; Daddi et al. 2007; Pannella et al. 2009; El- 
baz et al. 2011). Alternatively, galaxy spectral energy distri- 
butions (SEDs) can be fitted to theoretical star formation 
histories (SFHs) to derive SFRs (e.g. Shapley et al. 2005; 
Law et al. 2012a). Since dust is included as a free parameter 
in these SFHs, intrinsic (i.e. dust-corrected) SFRs can be 
derived self-consistently using this method. The derivation 
of reliable SFRs ideally requires rest-frame UV photometry 
(as is the case here), since the leverage in the SFR and the 
dust extinction comes largely from these wavelengths. 

In this paper, we calculate galaxy SFRs via SED fitting. 
The WFC3/ACS photometry of each individual galaxy is 
compared to a large library of synthetic photometry, con- 
structed using constant SFHs, each described by a stellar 
mass (M), age (T), metallicity ( Z ) and internal extinction 
( Eb-v )• We vary T between 0.05 Gyrs and the look-back 




Figure 2. TOP: SFR as a function of galaxy stellar mass in 
our galaxy sample. The main sequence of star-forming galax- 
ies at z ~ 2 is shown using the solid line (the observed scat- 
ter is indicated using the dotted lines). Galaxy morphologies are 
shown colour-coded (black = spheroids, blue = non-interacting 
late-types, red = major mergers). BOTTOM: The distribution of 
derived extinction values for our galaxies. The median value is 
shown by the vertical red line and the median from the recent lit- 
erature (see Tresse et al. 2007 and Cucciati et al. 2012) is shown 
by the blue dotted line. 


time to z = 20 in the rest-frame of the galaxy, Z between 
0.1 Z© and 2.5 Z© and Bb-v between 0 and 1 mag. Syn- 
thetic magnitudes are generated by folding the model SFHs 
with the stellar models of Bruzual & Chariot (2003), with 
dust attenuation applied following Calzetti et al. (2000). The 
likelihood of each model, exp(— x 2 /2), is calculated using 
the value of x 2 i computed in the standard way. Estimates 
for parameters such as stellar mass, internal extinction and 
SFR are derived by marginalising each parameter from the 
joint probability distribution, to extract its one-dimensional 
probability density function (PDF). We use the median of 
this PDF as the best estimate of the parameter in question, 
with the 25 and 75 percentile values (which enclose 50% 
of the probability) yielding an associated uncertainty. Since 
dust is explicitly taken into account in this process, we de- 
rive intrinsic SFRs, free of internal reddening, directly from 
the SED fitting process. The derived internal extinctions, 
SFRs and stellar masses are uncertain by ~0.1 mag, ~0.I 
dex and ~0.2 dex respectively. 

In the top panel of Figure 2, we plot our derived SFRs 
vs. galaxy stellar mass. The SFR values for our star-forming 
galaxies are consistent with the Star formation main se- 
quence at these epochs defined by the recent literature. 
Galaxies that lie below this sequence are typically spheroids, 
which are relatively passively-evolving systems. Recall that, 
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unlike studies that specifically target star-forming systems, 
the mass-complete sample employed here is not biased 
against galaxies with low star formation rates. In the bot- 
tom panel of Figure 2, we present the distribution of derived 
internal Eb-v values for our galaxies. The spread in our val- 
ues (0 < Eb-v < 0.5 mag) agrees well with that found by 
other studies (e.g. Law et al. 2012b) and the median of our 
distribution ( Eb-v ~ 0.25 mag i.e. Afuv ~ 2.2 assuming 
Calzetti et al. 2000) is in good agreement with the literature 
at z ~ 2 (see e.g. Tresse et al. (2007), Law et al. (2012b) 
and Cucciati et al. (2012, see their Figure 4)). 

4 THE MAJOR-MERGER CONTRIBUTION 

TO THE STAR FORMATION BUDGET 

We begin by exploring how star formation activity is ap- 
portioned in terms of galaxy morphology, by summing the 
derived SFRs of galaxies in each morphological class and 
considering the fractional contribution of these classes to the 
total star formation budget (Figure 3). We find that ~55% 
of the star formation activity takes place in non-interacting 
late-types, with 27% in major mergers and the rest (18%) in 
systems that have spheroidal morphology. It is worth noting 
that the proportion of star formation driven by morpholog- 
ically selected major mergers calculated here is higher than 
the corresponding value derived for starbursts by Rodighiero 
et al, (2011). As we noted in the introduction, this is due 
to the fact that many major mergers exhibit similar or only 
modestly-enhanced SFRs compared to normal star-forming 
galaxies, he on or close to the star-forming main sequence 
(see Figure 2 above) and are, therefore, not part of the more 
extreme starburst population. 

The predominance of non-interacting late-types in the 
total star formation budget indicates that major mergers 
are not the dominant mechanism driving star formation in 
massive galaxies at z ~ 2. Furthermore, as we noted in the 
introduction, systems undergoing major mergers continue to 
experience star formation via other processes (e.g. cold flows 
and minor mergers). Hence, 27% represents an upper limit 
to the major-merger contribution to the star formation bud- 
get. To improve our estimate, we consider the enhancement 
of star formation due to major merging, since this better 
represents the portion of the star formation activity that is 
directly attributable to this process. While measuring this 
enhancement is not possible in individual major mergers, we 
car. estimate a typical value for the population as whole by 
comparing the mean specific SFR in the major mergers to 
that in the non-interacting late-types. 

The ratio of the mean specific SFRs in these two mor- 
phological classes is ~2.2:1 (major mergers : non-interacting 
late-types), implying that, on average , around half the star 
formation in major-mergers are likely driven by other pro- 
cesses. This relatively modest enhancement in star formation 
activity due to major merging is consistent with the findings 
of recent theoretical work (e.g. Cen 2011) and also empiri- 
cal studies that do not find significant differences between 
the SFRs of galaxies that are morphologically disturbed and 
those that are not at this epoch (e.g. Kaviraj et al. 2012; 
Law et al. 2012a). Thus, if around half the star formation 
in major mergers is unrelated to the merger itself, then the 
major-merger contribution to the total star formation bud- 
get is likely to be as low as ~15% (i.e. 27% X 1. 2/2.2). 

Before we conclude this section, we briefly discuss the 
spheroid population in the context of the major mergers. We 
note first that the time inter\al spanned by our study (~0.3 


Morphology Fraction of SF budget 

Non-interacting late-types O.SS* 014 

Major merger 0.27^°®® 

Spheroids 0.18 ±OOB 
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( 27 *) 



Figure 3. TOP: The fraction of the total star formation budget 
in massive (M„ > 10 10 Mg) systems at z ~ 2 that is hosted by 
various morphological types. Major mergers account for less than a 
third (27%) of the total star formation budget, while non-interacting 
late-type galaxies host more than half of the star formation activity. 
Given that systems undergoing major mergers continue to experience 
star formation driven by other processes, 27% is an upper limit to 
the major merger contribution to the star formation budget. The 
actual contribution is likely to be ~15% of the total star formation 
budget (see text in Section 4 for details). BOTTOM: A pie-chart 
visualisation of the star formation budget apportioned in terms of 
galaxy morphology (Sph = Spheroids, LTG = Non-interacting late- 
type galaxies, Maj. Mgr — Major mergers). 

Gyr) is shorter than the effective timescales (0.5-2 Gyr) over 
which major mergers coalesce (see e.g. Lotz et al. 2008; New- 
man et al. 2012), so that the morphological classes do not 
overlap with each other. More importantly, however, Kaviraj 
et al. (2012) have used high-resolution cosmological simula- 
tions to demonstrate that spheroids at 1 < z < 3 that are 
remnants of recent major mergers (i.e. ones that coalesced 
within the last ~0.5 Gyr) will exhibit clear tidal features 
at the depth of the ERS images. This study has further 
demonstrated that many newborn spheroids in this redshift 
range do not carry such morphological disturbances, indicat- 
ing that a significant fraction of these systems are not built 
via major mergers (in agreement with the results of recent 
theoretical work, e.g. Dekel et al. 2009). Around 15% of the 
spheroids in our sample show morphological disturbances 
and these galaxies account for ~3% of the total star forma- 
tion budget. While the spheroid arid major merger classes 
do not overlap (as discussed above), it is clear that, even if 
we added the disturbed spheroids to the major merger por- 
tion of the star formation budget, our conclusions would re- 
main unchanged. Our analysis therefore indicates that major 
mergers contribute a relatively insignificant fraction (~ 15 %) 
of the total star formation budget in massive galaxies at 
zee 2 and are not the principal driver of cosmic star format 
tion at this epoch. 

5 SUMMARY 

We have explored the significance of major mergers in driv- 
ing star formation at high redshift, by quantifying the con- 
tribution of this process to the total star formation budget 
in a sample of 80 massive (M* > 10 10 Mg) galaxies at z re 2. 

We have found that ~55% of the total star formation activ- 
ity in massive galaxies at this epoch is hosted by late-type 
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galaxies that are not interacting with other systems, with 
~27% in major mergers, and the rest in spheroids. 

Since systems undergoing major mergers continue to 
experience star formation driven by other processes (e.g. 
cold flows and minor mergers), 27% is likely to be an up- 
per limit to the contribution of major mergers to the total 
star formation budget. To obtain a more accurate estimate, 
we have considered the typical enhancement of star format 
tion induced by a major merger, since this is the portion 
of star formation activity which is directly attributable to 
thi: process. We have estimated this enhancement using the 
ratio of the mean specific star formation rate in the major 
merger population to that in the non-interacting late-type 
galaxies. In agreement with recent observational and the- 
oretical work, we have found a relatively modest enhance- 
ment (around a factor of 2), which implies that, on average, 
roughly half the star formation activity in major mergers 
is unrelated to the merger itself. This reduces the contribu- 
tion of major mergers to the total star formation budget in 
massive galaxies at this epoch to ~15%. Our analysis there- 
fore indicates that the contribution of major mergers to the 
total star formation budget in massive galaxies at z ~ 2 is 
relatively insignificant and this process is not the principal 
driver of cosmic star formation at this epoch. 

Since our study is based on 'instantaneous’ star forma- 
tion rates, it provides only a snapshot of the star format 
tion budget in massive galaxies at z ~ 2. Thus, while our 
data cannot rule out a merger-dominated era in the very 
early Universe, if the merger contribution to stellar mass 
assembly does not evolve significantly into earlier look-back 
times (z > 2), then major mergers are unlikely to be sig- 
nificant contributors to the overall buildup of stellar mass 
in the Universe. In forthcoming papers we will use morpho- 
logical analyses of large datasets such as CANDELS - e.g. 
via projects such as Galaxy Zoo (Lintott et al. 2008), which 
uses 450,000+ members of the public to visually classify 
large survey datasets - to comprehensively study the con- 
tribution of major mergers to the star formation budget as 
a function of stellar mass, environment and redshift. 
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